CMB anisotropies seen by an off-center observer in a spherically symmetric 
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The current authors have previously shown that inhomogeneous, but spherically symmetric uni- 
verse models containing only matter can yield a very good fit to the SNIa data and the position 
of the first CMB peak. In this work we examine how far away from the center of inhomogeneity 
the observer can be located in these models and still fit the data well. Furthermore, we investigate 
whether such an off-center location can explain the observed alignment of the lowest multipoles of 
the CMB map. We find that the observer has to be located within a radius of ~ 15 Mpc from the 
center for the induced dipole to be less than that observed by the COBE satellite. But for such 
small displacements from the center, the induced quadru- and octopoles turn out to be insufficiently 
large to explain the alignment. 
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I. INTRODUCTION 

In a recent work [l[ , we studied spherically symmetric in- 
homogeneous universe models - the so-called Lemaitre- 
Tolman-Bondi (LTB) models. We found that for a cer- 
tain class of inhomogcncities, such models could easily ex- 
plain various cosmological observations without introduc- 
ing dark energy, most notably the luminosity distance- 
redshift relation of type IA supernovae and the position 
of the first peak in the CMB spectrum. The inhomo- 
geneities required are of the form of a spherically sym- 
metric underdense bubble in an otherwise flat and ho- 
mogeneous Einstcin-de Sitter universe, with the observer 
located at the center of the bubble. 

Unless the observer is positioned exactly at the center of 
the bubble, the distribution of matter, as seen by the ob- 
server, will be anisotropic. This will affect the observed 
microwave background and constrain the possible loca- 
tion of the observer, since the CMB dipole must be in 
agreement with observations Q. Note that in a homo- 
geneous universe model, this dipole is attributed to the 
peculiar velocity of the observer. However, as discussed 
in 3, in an LTB model there will be an additional con- 
tribution to the dipole from the anisotropy of space-time. 
Thus, the dipole seen by an off-center observer will be due 
to a combination of kinematic effects and the off-center 
location. 

The anisotropy will also induce higher multipoles in the 
CMB spectrum. Moffat Q proposes this mechanism as 



a possible explanation for the observed alignment of the 
CMB quadru- and octopole 0,0,0,1,11, 

since the direc- 
tion from the observer towards the center of the bubble 
singles out a "special" axis. 

In this work we will investigate these induced anisotropics 
in the CMB to establish how far from the center the ob- 
server can be located, and whether they can offer an ex- 
planation to the alignment of the lowest multipoles. We 
find that the observer has to be located within a sphere 
extending approximately 15 Mpc about the origin, in or- 
der for the induced dipole to remain within the observed 
range. However, within this small volume the induced 
quadru- and octopole turn out to be insufficiently large 
to explain the alignment. 

The paper is organized as follows: In Sect. [Til wc deduce 
the differential equations governing the path and rcdshift 
of photons. We then proceed to work out expressions for 
the induced temperature distribution and corresponding 
CMB multipoles in Sect. [Til] Next, in Sect. 1171 we solve 
these equations numerically to find the multipoles as a 
function of the position of the observer. This allows us 
to find LTB models which agree with both the observed 
dipole and observations of SNIa and the position of the 
first CMB peak, as described in [l[ . We present two such 
models in this section. Finally, in Sect. [V] we discuss and 
summarize our work. 
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II. THE GEODESIC EQUATION IN THE LTB y 
SPACE-TIME 



Th e sp herically symmetric LTB metric is given by fToL 

mm 



-dt 2 



[R'(r,t)]< 
1 + (3{r) 



dr 2 + R 2 (r,t)dn 2 , (1) 



where i?(r, t) is a position-dependent scale factor, and 
P(r) is related to the curvature. 

One of the great advantages of working in this space-time 
is that the Einstein equations can be solved exactly in 
the matter-dominated scenario. The function R(r, t) can 
then be written in terms of a conformal time 77, defined 
by (3 1/2 dt = Rdrj, as 

R = 




— (sinh 77 - rf) 



+Ro 



sinh 7] 




sinh rj 



(cosh rj — 1) 



(2) 



( (3) 



where a(r) > is an arbitrary function, and we have 
assumed j3(r) > 0. Furthermore, we have defined Rq = 
R(r, 0), with the initial time t = defined as the time of 
last scattering. 

Photons follow trajectories determined by the geodesic 
equation, 



d 2 x^ 
~dX r 



dx a dx v 
'~dX~dX 



0. 



(4) 



where T^ v is the Christoffel symbol, and A is a monoton- 
ically increasing (or decreasing) parameter defined along 
the path of the photons. 

Due to axial symmetry, the photon paths must be in- 
dependent of the azimuth angle (f>, which leaves three 
possible choices for free index \x. First, fi — t yields 



dX 2 



R'R' f dr 



dX 



1 + !3 
Next, jj, = r yields 

d 2 r (R^_ 
d\ 2 + \R r 2- 

2R! dr dt 
+ IFdXdX 

and finally, fi = 9 yields 



dQ 

RR[ dx 



0. 



(5) 



2/3 

_ RQ. + 0) 
R' 



= 0. 



f0_ 

d\ 2 



R' d9 dr R d9 dt 



R dX dX 



RdX dX 



(6) 



(7) 




FIG. 1: A photon hitting the observer at an angle £ 



which can be written as conservation of angular momen- 
tum J 



d\ 



R 



,d0 
dX 



= Tx J = - 



(8) 



In addition, the 4- velocity identity, u^u^ = for photons, 
leads to the constraint 



(it 



M2 



dr 



1 + \dX 



R 2 



(9) 



It is simplest to specify the initial conditions at the time 
to when the photon arrives at the observer's position, 
which is given by r = r$ and 8 = 0. The path of the 
photon is shown in Fig. [TJ It hits the observer at an 
angle £ relative to the z-axis. The spatial components of 
the unit vector along this axis are 



,• = ^2(1,0,0) 



(10) 



where the three components are in the r, 8 and <p direc- 
tion, respectively. 

The spatial direction u l is given by the tangent to the 
photon path at to, i.e. 



dX 
~dt 



dr dd dd>\ 1 , , „ , , „ 

dx>Tx>Tx J—fc'/tf.o). (ii) 



where the first factor ensures normalization, gijtfu 3 = 1, 
and we have introduced u = dt/dX and p = dr/dX. Note 
that we have chosen to let A decrease with time, since we 
will start integrating the equations at t = to and follow 
the photons backwards in time until recombination. 



The angle £ is given by the inner product of v l and u 

R' p 



cos£ = gijU l v J 



(12) 
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Since the parameterization of the photon path in terms 
of the affine parameter A is arbitrary, we can choose this 
such that A = when t — to and uq = u(X = 0) = — 1. 
Using Eqs. © and l|12p. these conditions translate into 
the following initial conditions 



Pa = — ^ — cos£. 
Jo = J = R sin £ . 



(13) 
(14) 



Furthermore, we need to determine the rcdshift of the 
incoming photons as a function of the direction. The 
photons follow a path given by r(A), t(X) and 8(X). The 
redshift can be represented by the change in time sepa- 
ration between adjacent photons. Due to the expansion 
of the universe, this separation changes as the photons 
propagate through space. The redshift is given by the 
relative change of the separation, i.e. z = (r r — r e )/T e , 
where the subscripts r and e refer to the receiver and 
emitter positions, respectively. 

Consider two photons emitted by a source with a time 
separation of r. Let the equation describing the time 
coordinate along the first geodesic be ti(A) = t(A). The 
time along the second geodesic is then given by £2 (A) = 
t(X) + t(A). Both photons must satisfy the 4- velocity 
identity For the first photon this reads as 



R'{r,tf 
1 + 



dr\ 2 , v2 fd9^ 2 



(15) 



while for the second photon, the 4-velocity identity be- 
comes 



d{t + T)Y R'{r,t + rf 



dX 



1 + /3 



(16) 



Expanding Eq. (|TB)) to first order in r and using Eq. (fTS")) , 

we arrive at the expression 



dt dr 
dXdX 



"(A) 



R'R' 
1 + 



RR 



(17) 



The redshift measured by the observer as a function of 
the period at the time of emission is defined as 



l + z(X e ) 



TjXr) 

r(X e ) 



(18) 



We differentiate this with respect to A e , which gives us 
the expression 



dz 
dX e 



1 dr{X e )T{X r ) 



r(A e ) dX e r(A e ) ^ ^ 

Next, using Eqs. (JTTJ) and (fTHJ) , and suppressing the sub- 
script e, we arrive at the equation 



dz 
dX 



-(1 + *) 



dA 

lit 



R'R' 



RR 



(20) 



This equation determines the change in redshift mea- 
sured by the observer along an infinitesimal distance dX. 
To find the redshift as a function of A for a photon hit- 
ting the observer today, we can sum up the infinitesimal 
contributions along the past light cone, 



cUn(l + z) 
dX 



R'R' n 



R 



(21) 



with the initial condition z(X — 0) = zq = 0. 



To summarize, we must solve the five first-order differen- 
tial equations for (t,r,8,p and z), with the correspond- 
ing initial conditions (to, ?"0i 0,po and zq), under the con- 
straints © and (p 



III. TEMPERATURE ANISOTROPIES 



We wish to examine how being situated away from the 
center of the LTB coordinate system affects the CMB 
temperature measured by the observer. Since space-time 
is no longer spherically symmetric around such an ob- 
server, we expect him to measure additional anisotropics 
in the temperature to those measured by an observer 
at the center. In this paper we concentrate on the ad- 
ditional anisotropies rising from the observer's location, 
i.e. we disregard any intrinsic anisotropics in the CMB 
temperature at the last-scattering surface. Thus, we as- 
sume the temperature at the last-scattering surface to be 
isotropic. Any anisotropics measured by observers today 
are therefore due to the propagation of photons through 
an anisotropic space-time. 

The temperature of the background radiation in a given 
direction is determined by measuring the intensity of inci- 
dent photons from this direction. Assuming the radiation 
to be black-body radiation, the intensity will be given by 
a Planck spectrum with a corresponding characteristic 
temperature. It can be shown that the radiation field 
preserves its black-body nature as it propagates freely 
through space under the influence of a gravitational field 
[l3| . At any later time, the spectrum will still remain a 
Planck spectrum, but with a different temperature. 

In our specific case, the CMB temperature seen today by 
an off-center observer is given by 



T(0 = 



r. 



(22) 



where T* is the temperature at the last-scattering sur- 
face, and £ is the angle defined in Fig. [TJ The average 
temperature T measured by the observer is then 



-in 



1 + Z(£) 



(23) 
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According to measurements made by the COBE satellite 
[3], this temperature is T = 2.725. We can now use 
Eqs. (j22|) and (|23|) to define an average redshift to the 
last-scattering surface: 



1 + 2* 



f 



1 



sin£ 



(24) 



The relative temperature variation measured by the ob- 
server today will then be 

e(0-^ = ^t^ = ^#. (25) 

K) T T 1 + ^(0 

It is often more interesting to consider contributions at 
different angular scales rather than the total anisotropy 
itself. Such an analysis can be performed by decomposing 
the temperature field in spherical harmonics Yi m : 



e (0 =X]<w^ 



(26) 



i. 



where the amplitudes in the expansion are recovered as 
o Jo 



a,, 



(27) 



These measure the level of anisotropy at different angu- 
lar scales, with larger I values corresponding to smaller 
scales. Since the relative temperature field does not de- 
pend on the azimuth angle </>, all the ai m will vanish, 
except those with m = 0. 

The observed dipole in the CMB is of the order |aio| ~ 
10 -3 . This will put a natural constraint on how far away 
from the origin the observer can be located, since a far- 
ther off-center position usually means a larger dipole. 



IV. NUMERICAL RESULTS 

Following [l[ we will consider LTB models where the in- 
homogeneity is an underdense bubble with a flat and 
homogeneous space-time outside. We will consider two 
specific models of this type corresponding to two differ- 
ent choices of the functions a(r) and j3(r) in Eqs. @ and 
d3]). We will refer to these as model I and model II. We 
parametrize these function in the same way as in [![, i.e. 
we write 



a(r) 
[3(r) 



= Htr 2 



'0 



11 r 

an — Aq ( tanh — . 

1 2 2 2Ar 



B - A3 ( - - - tanh L— — 
1 2 2 2Ar 



(28) 
(29) 



which corresponds to a smooth interpolation between 
two homogeneous regions (i.e. a spherical bubble in an 
otherwise homogeneous universe). The parameter Hq = 




■Model I 
-Model II 



1000 2000 3000 4000 5000 6000 



Physical distance [Mpc] 



FIG. 2: The matter density today as a function of physical 
distance for the two models considered. 



100 /i out km s 1 Mpc 1 is the Hubble constant of the outer 
homogeneous region today, while an and Bo = 1 — an 
are the relative densities of matter and curvature in this 
region. Further, Aa = —A/3 determines the difference 
in matter density between the regions, while rn and Ar 
specify the position and width of the transition. 

In our previous work [l[, we assumed that the observer 
was positioned at the center of the bubble, and found 
a model that gave a good agreement with the Hubble 
diagram of observed SNIa and the position of the first 
CMB peak. Model I is identical to this model, while 
model II is slightly different. The matter distribution to- 
day for these models is plotted in Fig. [51 where we have 
used the generalized matter density defined in pj. Var- 
ious other properties of these are listed in Tab. |U One 
notable difference between these two models is that the 
transition from the underdensity to the homogeneous re- 
gion is much sharper in the second model. Note that the 
physical values given in the table are found assuming that 
the observer is placed at the center. The shift parameter 
iS is defined in [lj , and is simply the shift of the first peak 
in the CMB power spectrum relative to the concordance 
ACDM model. As we can see, both models yield a very 
good fit to both SNIa and the first CMB peak. 

Figs. [3] and |4] show the solutions of the geodesic equations 
for two off-center observers in the two models, with the 
observers located 20 Mpc and 200 Mpc from the center. 
The blue lines show the photon paths in the metric (r, 9) 
space, while the red circles are the positions of incoming 
photons at evenly spaced points in (cosmic) time. The 
distortion of light paths is clearly visible on the bottom 
plot of Fig. 21 in which the observer is placed relatively 
far from the center of the inhomogeneity and the density 
gradient is large at the transition. 

In our calculations, recombination is assumed to occur at 
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Description Symbol Model I Model II 



Density contrast parameter 


An 


0.90 


0.78 


Transition point [Gpc] 


''(1 


1.34 


1.68 


Transition width 


Ad/d 


0.40 


0.03 


Fit to supemovac 


v 2 

XSN 


176.2 


177.8 


Position of first CMB peak 


S 


1.006 


1.002 


Age of the universe [Gyr] 


to 


12.8 


12.7 


Relative density at the center 




0.20 


0.25 


Relative density outside undcrdensity 


n m . out 


1.00 


1.00 


Hubble parameter at the center 




0.65 


0.63 


Hubble parameter outside undcrdensity 


h out 


0.51 


0.51 



TABLE I: The parameters and features of two inhomogeneous 
models, with the observer placed at the center. Note that do 
and Ad are physical distances corresponding to ro and Ar 




r cose 



FIG. 3: Model I: Photon paths converging on the observer 
located 20Mpc (top) and 200 Mpc (bottom) from the center. 
The red cross marks the center of the underdensity, while the 
red circles show positions that are evenly spaced in cosmic 
time with a separation of 1 Gyr 



t = and today (to) is defined to be the time when the 
redshift of photons emitted at t — reaches ~ 1100. 
(The exact value of depends on the matter density 
outside the bubble, and a fitting formula is given in [HI). 

As discussed in the previous section, an off-center ob- 
server will measure a temperature anisotropy due to the 
non-symmetric paths traversed by CMB photons in dif- 
ferent direction in the sky. Using Eqs. (|2TJ)) and (|2"T)) . 




r cose 



FIG. 4: Model II: Photon paths converging on the observer 
located 20 Mpc (top) and 200 Mpc (bottom) from the center. 
The red cross marks the center of the underdensity, while the 
red circles show positions that are evenly spaced in cosmic 
time with a separation of 1 Gyr 



we can now calculate the temperature multipoles seen by 
such an observer. As an example, a plot of the multi- 
poles can be seen in Fig. [5] for an observer who is located 
200 Mpc from the center in model I. 

In Fig. [51 the coefficients a; for the dipole (I = 1), 
quadrupole (I = 2) and octopole (I = 3) are plotted 
as functions of the observer's position in model I. The 
most striking feature of these plots is that the quadru- 
and octopoles are very small compared to the dipole. If 
we assume that the induced dipole must be smaller than 
10 -3 , the induced quadrupole is less than 10 -7 while the 
induced octopole is smaller than 10 -9 . 

In Fig. the a;o's of model II are plotted as functions 
of the observer's position. It is evident that the behavior 
is almost indistinguishable from that of Model I, except 
for the largest distances where the transition from un- 
derdensity to flat space is starting to show in the first 
model. 



6 







-2.5 -2 -1.5 -1 -0.5 0.5 1 



FIG. 5: Temperature anisotropy in galactic coordinates seen by an observer 200 Mpc from the center of Model I, oriented so 
that the induced dipole coincides with the direction of the dipole seen by the COBE satellite [lj. The top plot to the left 
shows the temperature map, which is completely dominated by the dipole. Then follow plots with the dipole, quadrupole and 
octopole removed, successively 



50 100 150 200 

Distance from origin [Mpc] 



Distance from origin [Mpc] 



Distance from origin [Mpc] 



FIG. 6: The aio as a function of the observer's position, in Model I. The dotted lines are linear, quadratic and cubic fits, 
respectively. 



10 100 150 200 

Distance from origin [Mpc] 



Distance from origin [Mpc] 



Distance from origin [Mpc] 



FIG. 7: The aio as a function of the observer's position, in Model II. The dotted lines are linear, quadratic and cubic fits, 
respectively. 
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DISCUSSION 



The main purpose of this paper has been to determine the 
maximum displacement of the observer from the origin of 
the undcrdcnsity, for which the induced CMB dipole re- 
mains in agreement with the results observed by COBE 
0. Of course, one could in principle introduce an ad- 
ditional peculiar velocity towards the center of the un- 
derdensity to compensate for a too large induced dipole, 
but such a coincidence would be very difficult to justify. 
Therefore, we must require the induced aio to be of order 
10 -3 or less, which from the plots in Figs. [H] and [7] can 
be translated to 



Lobs 



< 15Mpc. 



(30) 



where d is the physical distance. When compared to 
the size of the underdensity, which according to Fig. [2] 
is around 1 500 Mpc, this means that if we are placed 
at a random position inside the bubble, there is roughly 
a chance of 1 to 10 6 that we end up inside the region 
allowed by Eq. (f3"0"|) . This is a rather strong violation of 
the Copernican principle, which states that we are not 
situated at a special place in the universe. On the other 
hand, a 10~ 6 probability is still much better than the 
infinitely improbable case of the observer being exactly 
at the center of the underdensity. Note that the size of 
the undcrdcnsity is dictated by the fit to the CMB and 
SNIa data. We have not been able to find smaller bubbles 
that fit these data as well as the models considered here. 

From Figs. [S] and [7| we see that the induced multipolcs 
become larger the farther away from the origin the ob- 
server is located, as we would expect. Thus, the largest 
possible quadru- and octopolcs with a dipole compati- 
ble with COBE measurements are those for an observer 
about 15 Mpc from the origin. However, at this rela- 
tively small distance, the values for these are of the order 
10 -7 for the quadrupole, and 10 -9 for the octopole. It 
is therefore clear that the induced quadru- and octopole 
cannot explain the observed alignment of the low-1 multi- 
poles in the CMB, since their contributions are negligible 
compared to the observed anisotropics (which arc of or- 
der 10 -5 ). Furthermore, any off-center placement must 
necessarily result in axial symmetric contributions to the 
CMB spectrum. Even if such contributions were of the 
correct order, Rakic et al. [l6[ show that they are very 
unlikely to explain the alignment. 

The smallness of the induced multipolcs can be under- 
stood from a simplified Newtonian picture. Compared 
to the homogeneous case with a spatially constant Hub- 
ble parameter h ou t, the observer at d b s has a "peculiar 
velocity" of roughly 



lout 



Vp hit 

/ ~ ~c ~ 3000Mpc 



d , 



bs 



(31) 



with respect to the origin. In such a picture, the tem- 
perature anisotropics measured by the observer are at- 



tributed to a Doppler shift of the CMB photons due to 
his motion. The change in frequency can then be written 
as [13 



v e 1 — cos £ 



(32) 



where v a is the frequency measured by the observer, 
and u e is the frequency relative to a stationary back- 
ground. The temperature shift associated with this fre- 
quency shift of the CMB photons is 



T y/T^W 



(33) 



The average background temperature measured by the 
observer will then be 



Thus, the temperature anisotropy can be written as 
AT(0 2 f3 



0(0 



a fi±|) 1-/3 cos € 



1 



(35) 



Using Eq. (|27|) . the multipolcs can now be calculated in 
terms of the velocity of the observer. We find that the 
leading contribution to a;o is P l , 



am ~ 



obs 



(36) 



Thus, using this simplified Newtonian picture, we expect 
the dipole to scale linearly, the quadrupole quadratically, 
and the octopole cubically with the observers position. 
In Figs. [6] and [7] we have plotted the real dependence 
of these multipoles on the observers position along with 
a best-fit d l obs dependence. As we see, the Newtonian 
picture gives a very good description. The expressions 
for the three lowest multipoles to the lowest order in j3 is 



aio 
a 3 o 




obs 



d 2 

3000 Mpc I u obs ' 



16tt ( hin-hout \ ]3 
175 V 3000 Mpc J u obs ■ 



(37) 
(38) 
(39) 



Numerically, this approximation yields aio = 1-4 x 10~ 3 , 
a-2o = 5.2 x 10~ 7 and 030 = 1.2 X 10" 9 for an observer at 
dobs = 15 Mpc in Model I, whereas the exact values are 
am = 1-4 x 10~ 3 , a 20 = 2.3 x 10~ 7 and a 30 = 1.3 x 10~ 9 
respectively 

Eqs. ([3"T|) - (f3"9")) imply that it is impossible to obtain suffi- 
ciently large values for the quadru- and octopole as long 
as the dipole is within the limits set by the COBE data. 
Note, however, that Tomita [l8[ has previously found 
relatively large values for the quadrupole in more sim- 
plified bubble models (where two homogeneous regions 
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are separated by a massive comoving shell). It is unclear 
to us why the Newtonian approach fails for his models. 
We have attempted to reproduce his results by intro- 
ducing very narrow transition regions in our continuous 
model, but the results we get for the multipoles are of the 
same order as those quoted above. We therefore expect 
Eq. (|3"fj)) to be roughly correct for all models of our type. 

In our analysis so far we have only considered contri- 
butions to the multipoles from the off-center placement. 
There will of course be additional contributions from var- 
ious sources such as the intrinsic primordial temperature 
anisotropies, the integrated Sachs- Wolfe (ISW) effect [l9| 
and a non- vanishing peculiar velocity of the observer. We 
have seen that when the dipole is constrained by data, 
the quadru- and octopoles due to the off-center placement 
are considerably weaker than those observed in the CMB. 
A possible way to obtain stronger quadru- and dipoles is 
to place the observer farther away from the center, while 
allowing one or more of the effects mentioned above to 
cancel out the excessive contribution to the dipole. 

However, concerning the first two effects, it is clear that 
neither of these can achieve such cancellation. Although 
there is no way of measuring directly the intrinsic dipole, 
it is reasonable to assume that it is of the same order as 
the neighboring multipoles, which are of order 10 -5 Sim- 
ilarly, we expect the contribution to the dipole from the 
ISW effect to be of the same order as for the quadru- and 
octopolc. Therefore, it is very unlikely that these effects 
are responsible for a chance cancellation of an excessive 
contribution to the dipole from the off-center placement. 



A non-vanishing peculiar velocity can reduce the dipole 
to any desired value as long as the velocity is chosen 
large enough. However, multipoles due to such motion 
will have a hierarchical scaling similar to that which we 
showed in the Newtonian case. Thus, even if we manage 
to obtain values for the dipole and quadrupole of the 
correct order, the octopole would still be too weak. From 
this we can conclude that even when combined with other 
effects, the off-center placement cannot provide sufficient 
power to both the quadru- and octopole. 

In summary, LTB models like the ones listed in Table U 
are not ruled out on the basis of these results, but they 
do require a violation of the Copcrnican principle, since 
the observer would have to be located at a very special 
place. The volume within which the observer can be lo- 
cated is severely constrained by the size of the dipole 
induced by an off-center placement of the observer. As a 
consequence of this, the quadru- and octopole turn out 
to have insufficient power to explain the observed align- 
ment. However, the LTB models remain an exotic alter- 
native to dark energy as an explanation of the apparent 
accelerated expansion of the universe. 
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